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The time development of the reduced density matrix for a quantum oscillator damped by coupling 
it to an ohmic environment is calculated via an identity of the Debye- Waller form. Results obtained 
some years ago by Hakim and the author in the free-particle limit Q are thus recovered. The 
evolution of a free particle in a prepared initial state is examined, and a previously published 
exchange 0,0 is illuminated with figures showing no decoherence without dissipation. 
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I. INTRODUCTION 

In a recent paper in this journal, van Kampen0| has re- 
examined dissipation and noise in a quantum oscillator, 
treating it as a sub-system coupled to an environment. 
In working out his model, he has introduced some nice 
methods. They are slightly simplified and modified in 
this paper to revisit the closely related problem of quan- 
tum coherence and decoherence. 

It is worth emphasizing that the harmonic oscillator 
is particularly simple, so that the analysis given here is 
not generalizable to more interesting systems. However, 
it may make up in explicitness what it lacks in general- 
ity. Indeed, the ultimate aim of this research is to further 
clarify the "positivity problem" in time dependent quan- 
tum statistical equations. The present work is a first 
step in that direction, in the hope that deriving largely 
known results in a simple way will clear the path. 

This simplicity may offer amusement if not instruction 
to my long-time friends and colleagues Jim Langer and 
Pierre Hohcnberg, despite the burdens of their high of- 
fices. It is a pleasure and honor to dedicate the paper to 
them. 



II. MODEL AND PRELIMINARIES 

The model is described by the Hamiltonian 

h = \[p 2 + n 2 Q 2 } + i z k {Pi + "tlQk + KQoM)] 2 }, 

(i) 

where the oscillator labeled will be called "the sub- 
system" and the others "the environment." The quantum 
mechanical (and classical) equations of motion obtained 
from Eq. (1) are 

Qo + [nl + Y^(a k /uj k ) 2 ]Q + ^2a k Q k =0 



} k + tolQk + a k Q = 0. 



(2) 



Fourier transforming Eq.(2) and eliminating Q k , one ob- 
tains 



g 1 (lu + iO + )Qo — 0, where 



k 



(4) 



z is a complex variable, and the i0 + in Eq.(3) introduces 
the causal boundary condition. As pointed out in ref.0, 
where it is called G, g~ x {z) has zeros on the real z axis 
corresponding to the normal mode frequencies, <jj v , of the 
coupled system of oscillators, and Qq(uj v ) is the ampli- 
tude of the sub-system oscillator in the vth mode. 
Now, ohmic dissipation |7j requires 



(5) 



for lu less than an upper cut-off ia c . Substituting this 
form into Eq. (4) yields, for uj c w, 



- 1 (u + iO + ) = lu 2 



n 2 

"0 



-ni 



2r) 

T JO 



dcu- 



(w + i0+) 2 - u 2 



(6) 



demonstrating very explicitly that Eqs. (1) and (5) do 
indeed construct a linearly dissipative environment, with 
damping constant r\, without changing the system fre- 
quency fin. 

In ref. |j] , it is noted that the orthogonal normal mode 
transformation matrix X defined by 



Q k — ^ X kv q v , 



(7) 



where the q v $> are the normal co-ordinates, is obtainable 
from the Green function given in Eq. (4). The normal- 
izations 



Wi = i 



qI+J2qI = q 2 [i 



k 



and 

2 it, ,2 



a, 2 ) 2 ]=l, (8) 



show that the amplitudes corresponding to the mode v, 
and thus the matrix elements of the transformation, are 
given by 



(3) 



2 



Xkv — 



1 1 dff~ ] 
2u>„ du v 

Oik 



and 



Xq v . 



(9) 



One also learns from ref.|^| that Eq.(9) may be used to 
deftly perform sums over normal modes. The complex 
function g{z) has poles only on the real axis and thus the 
spectral representation 



g[z) 



dbJ s(ll>) 

2ir z — uj 



(10) 



Using the explicit form for Xq v given in the second line 
of Eq.(9) one obtains a formula that will be useful later 
in this work: 



E 



/dz 
—zg(z)F(z) 



= 2 £° ^ws(u,)F(u,). (11) 

The contour surrounds the real axis, where the function 
F has been assumed to be regular and zero for u> < 0; it 
has been evaluated using Eq.(10). 

Note also from Eq.(6) and its complex conjugate that 
for an ohmic environment 



s(uj) 



2u)T) 



(UJ 2 - 05) 2 + LU 2 7] 



(12) 



III. TIME EVOLUTION OF THE REDUCED 
DENSITY MATRIX 

To carry out the program of this section, initial condi- 
tions must be specified. I shall use those of and 0, 
which are a special case of ones originated, to the best of 
my knowledge, in ref. 1]. Assume that at t = complete 
thermal equilibrium is disturbed by a real "aperture func- 
tion" a(Qo). The entire system is then allowed to evolve 
to time t, and projected on to position states of the sub- 
system. The resulting reduced density matrix is given by 
(h = 1) 



Tr{\QZ f )(Q' of \e- 



76 f ,t) (13) 

:m a(Q )p th (H)a(Q Q )e +im }. 



In this equation, the primed quantities are ordinary num- 
bers, the unprimed ones operators, Tr indicates a trace 
over all states of H, and pth = exp{— /3i7}/Trexp{— f3H} 
with (3 the reciprocal temperature. Fourier transform the 
aperture function 



&(Qo) — / daa(a)e 



-iQaa 



(14) 



and express the projection operator in Eq.(13) as 

\Qof){Q'of\ = J dudvf(u,v)e iPoU e lQoV with 

f(u,v) = ±- e -^f v 6(Q' of -Q: f -u) 7 (15) 

proved by taking matrix elements of both sides, to see 
that 



p(Q'o f , Q'o f , i) — / dadbdudva(b)a(a)f(u,v)T, 



with T 



JQab iP (t)u iQ {t)v iQ a\ 



(16) 



Here the brackets mean an average with respect to pth, 
and operators with a time argument are in the Heisen- 
berg picture. A somewhat simpler average is done fairly 
heroically in ref. 3 using properties of Laguerre polyno- 
mials. However, since the sub-system and environment 
are all harmonic, a low-brow method is available. A 
single simple harmonic oscillator obeys the well known 
Debye- Waller identity for thermal averages, 



w\ = c -i c2 (<7 2 ). 



(e tqc ) 



(17) 



where q is the position operator and c a number. This is 
reviewed in the Appendix, where it is also shown that a 
straightforward generalization yields 



lnT= - i[( a + b) 2 +v 2 }{Ql)-±u 2 {P§)-bu{Q P Q (t)) 

- bv(Q Qo(t)) - uv(PoQo) - ua(P (t)Q ) 

- av(Q {t)Qo). (18) 

Since the co-ordinate Qo and the momentum Pq are lin- 
early related to the normal mode q u s and p u s via the 
known Ao^s the correlators in Eq.(18) are readily calcu- 
lable, thereby formally completing the task of this sec- 
tion. 



IV. AN EXAMPLE 

To illustrate the usefulness of these methods, consider 
(Qo), one of the averages occurring in Eq.(18). From 



Eq.(7) 



= E X ol(^cothi/3^), (19) 



because the us are independent oscillators, and (q v ) = 0. 
The sum can now be transformed using Eqs.(ll) and 
(12). It is easy to do analytically at zero temperature 
(j3 = 00). Define the real part of the damped oscillator 



3 



frequency via f7 2 = O 2 , — ry 2 /4 and factorize the denomi- 
nator in Eq. (12) to obtain 
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(20) 



The last is a known answer for the dissipation-induced 
squeezing by an ohmic bath. 0, @ This effect is at the 
root of the effect of dampin g in reducing the rate of es- 
cape from a metastable well.|ll| 
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FIG. 1: Time evolution at low temperature, lth/d = 1. [The 
symbols are defined in the text.] 



V. TIME DEPENDENCE OF THE PROBABILITY 

The position-space probability for the sub-system 
P(x, t) is obtained by setting Q' j — Q'^ — x in Eq.(16), 
whereupon Eq.(15) requires that the variable u =4- in 
Eq.(18). The time dependence can then be completely 
described by 



([Qo(t) - Qo(0)]Qo(0)) = -C(t) + iA(t), 



(21) 



this being the notation used in ref . Q in the free-particle 
limit, with the sign as corrected in footnote 0. One can 
now do the (gaussian) integrals in Eq.(16) to obtain, with 
x ' = Q'oi ~ Q'iii X' = 2 {Q'oi + Qoi)' 



P(x,t) = J dQ' oi dQ'^a(Q' 0i )a(Q'^)J(x,x',X',t), 

11 x' 
J = 7tt — exp \i — — — (x — X')] 



x exp — [ 



x' C(t) 1 .x'C{t) 



AA>{t) 2{Ql 



2A{t) 



.(22) 



In the free-particle limit flo — > 0, (Qq) — > oo and the 
normalized gaussian in X' must be replaced, to preserve 
dimensions and normalization, by where L is the 

size of the system, yielding 



J 



exp \i — -rr^-(x — X 



„ x' 2 C(t) 
4A 2 {t) 



4irA(t)L r 2A{t) 



(23) 



Using the method of Section II, known results for a free 
particle follow from Eq. (21): 



(I 



lu(uj 2 + i] 2 ) 2r\ 



r\ f Uc , , f3u (1 - cosujt) 
C{t) = - / dweoth-— ^—5 



(24) 



Wv » i) 



i 



l-e-^-rjt). 
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FIG. 2: Time evolution at moderate temperature, hh/d 



VI DECOHERENCE 

Eq.(22) — which is not restricted to ohmic dissipation — 
is to my knowledge new, and explicit enough to allow 
a general study of time development in this dissipative 
quantum system. Previous attempts Il2| in which I 
have been involved have for technical reasons been re- 
stricted to high temperature, a limitation which it should 
here be possible to avoid. Work, in progress, in this direc- 
tion would seem to be justified by experimental [l3| and 
theoretical|l4j interest in quantum information storage. 

The free particle limit, Eqs. (23,24), has been derived 
by many different methods — perhaps none as straightfor- 
ward as the one here given. At finite 7], these equations 
(uncontroversially) display decoherence. This is partic- 
ularly well demonstrated in an example introduced in 
ref [Tol]. Here the aperture function is taken to be a sum 
of two Gaussians, each of width a and separated by a 
distance d. The probability given by Eqs. (23) and (24) 
can then be written as a sum of the probabilities from 
each slit alone (sum term) and an interference contribu- 
tion depending sinusoidally on a time dependent phase 
(interference term.) There is an unresolved controversy 
0, Q m the published literature about what is meant 
by decoherence in this completely well defined problem. 
No one would doubt that the amplitude of the inter- 
ference term is a measure of coherence. In ref. and 
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FIG. 3: Time evolution at high temperature, hh/d —1/5 

other publications^^] an "attenuation coefficient" is in- 
troduced which is equivalent to the amplitude of the in- 
terference term divided by the sum term, evaluated at 
the mid-point between the slits. This quantity is by con- 
struction unity at t = 0. It decreases rapidly with time. 
At high temperatures, in a sense to be made clear below, 
it drops to zero, even when the environmental coupling 
is eliminated. This is interpreted in refs.|0,^^| as "De- 
coherence without Dissipation." 

In a Comment Q it is suggested that the measure 
of decoherence used in 0] does not distinguish be- 
tween the loss of phase information and the spreading 
of wave packets on time scales less than the mixing time 
tmix = 2ad(m/K). Here a is the width of each slit, d the 
spacing between them, and I have reintroduced Planck's 
constant Ti and the particle mass m to make the dimen- 
sions transparent. 

This has been vigorously rebutted in a Reply Q- 
Rather than repeat these arguments, I refer the in- 
terested reader to them. However, since many readers 
may be intrigued by the idea of decoherence without dis- 
sipation, I close this paper with 4 figures which show 
that there is no evidence for any such thing in the un- 
controversial Eqs. (23,24). Since the disagreement occurs 
in the limit of no dissipation, consider this case at vari- 
ous temperatures — given by the ratio of the thermal de 
Broglie wavelength l t h = V '/3 /2(?l/ \pm) to d. The fig- 
ures, in which a/d — 0.05, show that at low temper- 
atures, hh ~ d there is coherence without decoherence. 
At higher temperatures, there is no coherence at all, even 
on the short time scale 

cP" I (J 

tflo = x/8/3m— = 2t mix ■ — ■ — (25) 
d d d 

introduced in and thus nothing to decohere. In 
this limit, coherence is already destroyed by the Hakim- 
Ambegaokar initial condition. 

Added Note. After this paper was completed, I was made 
aware of in which normal co-ordinates are used to 
treat the problem of many non-interacting fermions cou- 
pled to a disspative environment in a harmonic oscillator. 




FIG. 4: Blow-up of Fig. 3 on the scale of tflo = 0.02 tmix, 
for the chosen parameters. 

APPENDIX 

A wonderfully short proof of the identity for thermal 
averages 

(exp[y^ djdj + Cjfflt]) = exp[i J]i((^ a » + c i a !) 2 )]> ( 26 ) 

i 

where a^, a\ are boson annihilation and creation opera- 
tors and the subscript i refers to independent harmonic 
oscillators, is given in a single mannerist, if not rococo, 
sentence by Mermin 16]. Note that, since all the oper- 
ators in Eq.(16) have c-number commutators, it can be 
put in the form of the left hand side of Eq.(26) for the 
subsystem oscillator labeled 0, using the Baker-Haussdorf 
identity for such operators: e A e B = e^ A+B ^e^ A,B]> . Now, 
express Pq and Qo in terms of the (normal) co-ordinates, 
of the independent oscillators v. The Debye- Waller iden- 
tity Eq. (26) then has in the exponent a sum over correla- 
tors for each v. Note that as in Eq. (19), a single sum can 
be replaced by a double sum, because {q v ) = (p„) = 0, 
yielding Eq.(18). 
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